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ABSTRACT 

^ ■ We report the discovery of a small cluster of massive stars embedded in a NIR 

^ ■ nebula in the direction of the IRAS 15411-5352 point source, which is related to 

^ ■ the alleged planetary nebula W16-185. The majority of the stars present large 

^ ■ NIR excess characteristic of young stellar objects and have bright counterparts in 

the Spitzer IRAC images; the most luminous star (IRSl) is the NIR counterpart 
of the IRAS source. We found very strong unresolved Br7 emission at the IRSl 
position and more diluted and extended emission across the continuum nebula. 
From the sizes and electron volume densities we concluded that they represent 



> 



^ ■ ultracompact and compact HII regions, respectively. Comparing the Br7 emis- 

■ sion with the 7 mm free-free emission, we estimated that the visual extinction 

^ . ranges between 14 and 20 mag. We found that only one star (IRSl) can provide 

O ■ the number of UV photons necessary to ionize the nebula. 
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1. Introduction 

W16-185 was classified as a planetary nebula by Wray (1966) but because of its prox- 
imity to the IRAS point source IRAS15411-5352, Acker et al. (1987) concluded that it is 
a compact HII region. Further investigation by Noumaru & Ogura (1993), through low- 
dispersion spectroscopy and narrow band CCD imaging in the Ha line and the adjacent con- 
tinuum, showed that the monochromatic Ha image has its centroid displaced by Aa = 11" 
and AS = —13" from the IRAS source although it lies within its error ellipse. Even so, since 
the object is located in the border between compact HII regions and PNe in the far infrared 
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color-color diagram (Pottasch et al. 1988), they suggested that it was better described as a 
low-excitation planetary nebula. 

Although there is no point source in the 5 GHz PMN catalog (Wright et al. 1994) in 
the direction of W16-185, the region can be considered part of the G326.7-I-0.6 radio source, 
associated with the RCW95 complex (Rodgers, Campbell & Whiteoak 1960). Goss &: Shaver 
(1970) mapped this region at 5GHz and found the maximum emission in the direction of the 
IRAS15408-5356 source and an unresolved extension in the direction of W16-185. Strong 
CS(2-1) and NH3(1,1) line emission were detected by Bronfman, Nyman & May (1996) and 
by Roman-Lopcs ct al. (2002) for both IRAS sources. The radial velocities of the CS and 
NH3 lines differ by less than 2 km s"^ between the two regions, suggesting that they are at 
the same heliocentric distance. 

Roman-Lopes & Abraham (2004) reported the existence of a rich cluster of OB stars 
centered at IRAS15408-5356. In this paper we report the detection of a newly formed 
cluster of stars embedded in a small NIR nebula, associated to the IRAS15411-5352 source. 
We observed the region in the J, H and K filters, as well as in the integrated Br7 hydrogen 
line. The continuum subtracted Br7 image shows two distinct regions, an unresolved ultra 
compact HII region (UCHII) and a more extended 10" compact HII region (CHII). Using 
the Itapetinga radiotclcscopc at 43 GHz, with 2'. 2 resolution, we were able to separate the 
radio source at the position of W16-185 from G326. 7-1-0.6. 

In Section 2 we describe the LNA and radio observations as well as the Spitzer/IRAC 
data. In Section 3 we identify the stellar cluster, the IRAS source and discuss the continuum 
and Br7 emission from the nebula. In Section 4 we summarize our results. 

2. Observations and data reduction 

2.1. LNA Near Infrared data 

The imaging observations were performed in three observing runs: the first two were 
performed in 2000 and 2003 May, using the Near Infrared Camera (CamlV) of the Labo- 
ratorio Nacional de Astrofisica (LNA), Brazil, equipped with a HAWAII 1024 x 1024 pixel 
HgCdTe array detector mounted on the 0.6 m Boiler & Chivens telescope. The images were 
made in the direction of IRAS point source 15411-5352, using the J and H broad band filters 
and the K narrow band filters: CI, centered in the continuum at 2.14//m and Br7, centered 
in the corresponding line. Those observing missions were part of a survey (Abraham et al. 
2003) aimed to detection of the ionizing stars of compact HII regions associated to IRAS 
point sources that are also strong CS emitters (Bronfman et al. 1996). 
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The resulting data showed the presence of several embedded NIR sources as well as 
intense Br7 emission in the direction of the IRAS source, however the K band image reso- 
lution obtained with the 0.6 m Boiler & Chivens telescope (plate scale of 0.47 arcsec pixel"^ 
and mean FWHM of about 2" ) limited our capability to study the stellar population in the 
direction of the nebula. In 2005 May we made new NIR imaging observations at LNA using 
the Near Infrared Camera attached to the 1.6 m Perkin- Elmer telescope. The plate scale 
of this new observation was 0.24 arcsec pixel""^ and the mean values of the point-spread 
function FWHM were 0.95", 1.0" and 1.1" at the J, H Bxid K { CI and Br-f filters) images 
, respectively. The total integration time for the J and H bands were 330 s and 210 s 
respectively, and 1190 s for the CI and Br^ filters. The completeness magnitude limits of 
17.5 (J), 16.6 {H) and 15.0 (CI) were derived from the points at which the number N(m) 
of detected sources with magnitude m deviate from a straight line in the log (A^) versus m 
diagram, as can be seen in Figure 1. The resulting frames were reduced using the IRAF 
reduction system and the Point Spread Function fitting photometry was performed using 
the DAOPHOT package (Stetson 1987). Details about the calibration and reduction proce- 
dures can be found in Roman- Lopes, Abraham & Lepine (2003). Despite the strong nebular 
contribution present in the nebulae region, which could introduce a bias to the photometry, 
we are confident in our PSF fits because the star subtraction process gave good results, with 
no over-subtraction against the nebulosity. 

To check the consistency of our photometry, we compared the LNA magnitudes with 
that of the 2MASS survey. First we evaluated the color correction term from 43 isolated 
common stars outside of the infrared nebula, with good H and iiTphotometry in the 2MASS 
survey (fiag A or B). As we can sec from figure 2 the color term correction for the LNA 
photometric system is very small. In addition, we made comparative magnitude diagrams 
for a sample of 101, 131 and 126 common stars in the J, H and K filters, respectively, which 
arc presented in Figure 3. We found a good linear relation between the two systems, with a 
slope of about 1 and a dispersion that increases with the magnitude. 

2.2. Radio data 

The 7 mm radio observations were made with the 13.6 m radome-enclosed Itapetinga 
radiotclcscope, which has a beam size of about 2'. 2 at this wavelength. The receiver is a X 
band room-tcmpcrature mixer with 1-GHz double side-band, giving a system temperature 
of about 700 K, and was operated in the total power mode. The observations consisted in 
a series of 20 s right ascension scans with 30' amphtude. A hnear baseline was subtracted 
from the data to eliminate the sky contribution. The atmospheric opacity was calculated 
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using the technique developed by Abraham & Kokubum (1992) and the flux density was 
calibrated by observation of the standard point source Virgo A. The quoted errors are the 
rms fluctuations of the flux densities, obtained by fitting a Gaussian function to the scans. 

2.3. Spitzer/IRAC data 

The Infrared Array Camera (IRAC) of the Spitzer Space Telescope^ operates at four 
near to mid-infrared bands centered at 3.6, 4.5, 5.8 and 8.0 /xm. The images used in this 
work are sections of the original frames of the SIRTF Galactic Plane Survey. The IRAG 
Post-Basic Calibrated Data (PBCD) of the region around IRAS1541 1-5352 were taken from 
the data archive of the Spitzer Science Center^. Photometry of the Spitzer sources were 
performed using DAGPHOT package within IRAF, with an on source aperture radius of 3 
pixels and a sky annulus extending from 3 to 7 pixels. The aperture corrections and zero 
magnitudes fluxes were applied following the procedure described in the Spitzer/IRAG user 
manual. 



3. Results &c Discussion 

In Figure 4 (upper left side) we show the ~ 3.9' x 3.5' combined false-color image, made 
from the observed J (green), H (green) and CI (red) images. We can see in the enlarged 
image (at the botton left side of that figure) several point sources embedded in a small and 
bright infrared nebula, suggesting the presence of a small cluster, which coincides with the 
position of the IRAS 15411-5352 error eUipse. Prom an inspection of the R band image of 
the region (taken from the digitized all sky survey - DSS^) we verified that the near infrared 
nebula lies just bellow the R image, implying that W16-185 is more extended than what is 
inferred from the optical data and that the extinction increases strongly to the south. 

As mentioned before, W16-185 can be considered part of the RGW95 complex; in the 
right part of Figure 4 we show the 5 GHz contour map of the radio source G326.6+0.6, 
associated to the RCW95 complex, obtained from Goss & Shaver (1970), in which the 
approximate position of the observed NIR images is delimited. Our 7 mm radio observations 
in the direction of IRAS15411-5352, with a 2.2' beamwidth, were able to separate the source 



^http:/ /irsa. ipac.caltech.edu/Missions/spitzer. html 
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at this position; in a 30' scan at constant declination, it appeared as an isolated Gaussian 
profile centered at the IRAS source position with HPW of 3.4'. The deconvolved angular size 
is 2.6' resulting in an integrated flux density of 5.2 ± 0.9 Jy. The 7 mm complete map of the 
RCW95 region will be published elseware (Barres de Almeida, Abraham & Roman-Lopes 
2005). 

3.1. The stelleir cluster 

In order to investigate quantitatively the star density in our images, we plotted in Figure 
5 the positions of all sources detected above 3 cr of the local sky in the LNA K band image, 
labeled as IRAS 1541 1-5352 panel. We can sec a high concentration of stars near the IRAS 
coordinate (the "nebula" field); yet a high density of sources can result from localized low 
extinction as discussed by Persi et al. (2000). On the other hand, young clusters embedded in 
their parental molecular cloud frequently show a high fraction of stars with excess emission 
at NIR wavelengths. To investigated the nature of the stellar population in the direction of 
IRAS1541 1-5352 field, we compared it with the 2MASS photometry of a 4.0' x 4.0' nearby 
field centered at a;=15h:44m:40s and 5=-54d:00m:00s, labeled as "control" in Figure 5, which 
probably consists of Galactic field stars. 

We analyzed the NIR colors of the two selected regions using the (J — H) versus (H — K) 
comparative diagrams shown in Figure 6. There, we also represented the position of the main 
sequence and the giant branch (Koornnccf 1983) as well as the reddening vectors; the location 
of 10, 20, 30 and 40 mag of visual extinction are also marked (Fitzpatrick 1999). 

In the "control" color-color diagram we can distinguish two stellar groups: the first 
is formed by foreground field stars and the other by cool giants or super-giants, probably 
located behind the molecular cloud. In the IRAS1541 1-5352 color-color diagram, besides 
the population found in the control region we can see other objects showing excess emission 
at 2.2 fim, probably due to the presence of warm circumstellar dust, characteristic of very 
young objects (Lada & Adams 1992). A significant part of this last group is located inside 
the region labeled "nebula" in Figure 5. 

In order to separate the cluster members detected only at the H and K bands from 
background field stars, we compared the number N{H — K) of objects with colors H — K 
with that found in the control population (normalized to the same area). The result is 
presented in Figure 7, where we can see that the foreground population is the same in both 
regions but there is an evident excess of sources with H—K > 0.5 mag in the IRAS1541 1-5352 
region. Using this result and the (J — H) x (H — K) diagram constructed from the sources 
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detected only in the nebula region (shown in Figure 8), we selected as member candidates 
all sources that present excess emission in the NIR, the objects with H — K > 0.5 that are 
located inside the nebula and do not have J — H colors of late-type stars and also all other 
sources with H — K > 1.5. In Table 1 we present the coordinates and photometry of all 
selected stars. 

Since the excess emission is generally attributed to the presence of disks or cocoons and 
background giant stars may contaminate the selected sample, we analyzed the recent images 
of the nebula made by the Spitzcr telescope, which are shown in Figure 9. The intensity of 
the extended emission that we sec in the Spitzer images increases with wavelength, implying 
that there is much more dust than what we can infer from the LNA images. Since Spitzer 
angular resolution is about 2 times smaller than that of our NIR observations, the sources 
inside the nebula region are not resolved, except for IRSl, which appears very bright and 
has the same coordinates in both surveys. The high extinction at the center of the cluster 
guarantees that the background contamination is small. 

To verify the nature of the LNA selected sources placed at the border of the nebula, we 
constructed their Spectral Enegy Distribution (SED) including the IRAC data and compared 
it with the SED of background red giants. In the late-type sample, the IRAC flux densities 

decrease with increasing wavelengths, while the opposite occurs with the cluster sources, 
suggesting the presence of material remaining from the star forming process. We illustrate 
this result in Figure 10, where we present the SED of IRS8, IRS19 and of late-type star. 

IRS7 and IRS16, appear as a single source in the Spitzer images although only IRS16 has 

excess emission in the J — H versus H — K diagram. IRS8, with H — K — 4.3, not detected 
at LNA J band and located at the west side of the NIR nebula appears as a bright IRAC 
source at 3.6 and 4.5/x. The same is valid for 1RS4, which does not show excess emission 
at 2.2yum and therefore can be classified as a BO-BO. 5V ZAMS. 1RS19 is one of the most 
interesting of the selected sources. In fact, not only was not detected at the LNA J band and 
has — X = 3.34, but as we can see from Figure 10, its flux density increases continuously 
from 2.2yum until 8.0y[/m. Finally there is another group of IRAC sources southeast of the 
NIR nebula, dehmited by the dashed box in Figure 9, from which only one was detected in 
the LNA image with magnitude K = 14.95±0.06 and coordinates (J2000) ci;=15h:45m:04.53s 
and (5=-54d:02m:39.6s. 

The luminosity of the objects in Table 1 can be derived from the J versus {J — H) 
color-magnitude diagram shown in Figure 10 where the locus of the main-sequence for class 
V stars at 2.4 kpc (Giveon ct al. 2002) is plotted and the position of the spectral types 
earlier than F8V are indicated by labels. The intrinsic colors were taken from Koornneef 
(1983) while the absolute J magnitudes were calculated from the absolute visual luminosity 
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for ZAMS taken from Hanson et al. (1997). The reddening vector for a ZAMS BO V star, 
taken from Fitzpatrick (1999), is shown by the dashed hne with the positions of visual 
extinction Ay = 10, 20 and 30 magnitudes indicated by filled circles. We also indicated 
the stars with and without excess in the color-color diagram by open and filled triangles, 
respectively. IRSl, the most luminous object in the region, presents large excess emission 
in the near infrared and is the best candidate to the near infrared counterpart of the IRAS 
source. 



3.2. The IRAS source 

In Figure 11 we show the flux calibrated contour diagram made from the CI band image 
centered in the IRAS 15411-5352 position with the error ellipse delimited by the broken line, 
and where we indicated by fllled circles all the cluster member candidates. Several sources fall 
inside the IRAS error ellipse; a more accurate coordinate for the IR source was obtained from 
the Midcourse Space Experiment - MSX^ point source catalog. The MSX surveyed the entire 
Galactic plane within | 6 |< 5° in four mid- infrared spectral bands centered at 8.28, 12.13, 
14.65 and 21.34 //m, with image resolution of 19 arcsec and a global absolute astrometric 
accuracy of about 1.9 arcsec (Price et al. 2001). We found one MSX source within the IRAS 
error ellipse, with coordinates a(J2000) = 15M4'"59.38^ 5(J2000) = -54'^02'^19.3^ which in 
fact are coincident with the star that we labeled IRS 1. 

In figure 12 we show the spectral energy distribution of the IRS 1 source between 1.25 
and 100 ^m, taken from the IRAS and MSX catalogues and from this work. We determined 
a lower limit for the spectral type of IRS 1 using the integrated near to far IR flux density, 
assuming that all the stellar luminosity is used to heat the dust that emits at infrared 
wavelengths, obtaining LboI > 7.7 x which corresponds to an 08.5V ZAMS star 

(Hanson et al. 1997). 

3.3. Br7 emission and the infrcired nebula 

In order to study the extended emission, we flrst scaled the By'j to the continuum image 
(CI filter) using the common bright field stars. Next we constructed the fiux calibrated 
contour diagram of the difference between the two images and, by measuring the area between 
contours, we obtained the net Br7 flux using the relation between magnitude and flux density 



^http: / / www.ipac.caltech.edu/ipac/msx/ msx.html 
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given in Koorneef (1983). We found very strong Br7 emission at the position of the IRAS 

source and more diluted and extended emission across the continuum nebula. The integrated 
Br7 fluxes are (9.1 ± 1.2) and (12.6 ± 1.6) in units of 10~^^ erg cm~^ s~^ for the compact 
and extended sources, respectively. 

In Figure 13 (left side), we show the Br7 contour diagram (white lines) together with 

the selected NIR sources (white stars) and the Ha contour diagram taken from Noumaru & 
Ogura (1993) (black lines), overlaying the CI LNA image. In the right side of the Figure we 
present the two dimensional Br7 intensity distribution, which clearly shows the compactness 
of the HII region associated to the IRAS1541 1-5352 source. In fact, the FWHM of this 
distribution is about 1.1" (~ 0.012 pc at the quoted distance), which coincides with the 
stellar PSF, showing that we still did not resolved the region. 

The same procedure used to integrate the Br7 emission was used to calculate the nebular 
fluxes at the J, H and K bands, resuhing in 0.10 ± 0.01 Jy 0.29 ± 0.05 Jy and 0.38 ± 0.04 
Jy, respectively. These values can be compared with the expected free-free emission from an 
optically thin plasma at wavelength A, given by: 



where E is the emission measure (cm~^), Q the solid angle of the source, T the temperature 
in Kelvin and gff{X,T) the dimensionless Gaunt factor, which is approximately 1 for near 
infrared and for radio can be calculated from (Lang 1978): 



A lower limit to the product Q,E can be obtained from the measured Br7 fluxes, not 
corrected by absorption, using the expression given by Watson et al. (1997): 



= 5.4 X W-^^T-'^-^gff{X,T)nEe-^'/^''^ Jy 



(1) 




(2) 



S'Br7 = 0.9/Z-I/Br7«B^7";~-^ ^^S^ Cm ^ S 



1 



(3) 



where 



a 



•Br7 



= 6.48 X 10~''T 



llrp-im 



(4) 



From these QE lower limits and assuming T = 7500 K, we computed lower limits for 
the free-free emission at the J, H and K bands, as well as for the radio emission at 7 mm. 
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which are presented in Table 2 together with the observed values. We can see that there are 
excess flux densities at the infrared bands, which decreases with the increase in frequency, as 
expected in the presence of wavelength dependent absorption. Since the absorption at the 
Br7 and narrow K bands are equal, the excess in this band can be attributed to thermal 
dust emission and/or scattered stellar light. 

From the measured sizes (0.12 pc for the extended source and less than 0.02 pc for 
the compact one) and the Br7 fluxes, we derived emission measures of 5.3 x 10^"* cm^"'' and 
3.8 X 10^^ cm^^, respectively. These results are lower limits because the measured Br7 fluxes 
were not corrected for extinction and, at least for the compact region, we only have an upper 
limit for its size. 

Using these values for the emission measure, we calculated the optical depth at radio 
wavelengths from the expression: 



r(A,r) = 2xlO-^3^Af^T0A^^^ (5) 

We obtained lower limits for r(7mm) of the order of 0.02 and 0.003 for the compact and 
extended regions, respectively. Notice that the value obtained for the optical depth of the 
extended region can increase by a factor of a few when Br7 absorption is included, but it will 
still remain much smaller than one. On the other hand, the true optical depth of the compact 
region can be much larger than 0.02, if the real size is smaller than our measured limit. To 
test this possibility, we derived the emission measure and the corresponding angular size 
under the assumption of T(7mm) = 1. We obtained £■ = 1.8 x 10^^ cm~^ and an angular 
size of 0.1", which is about eight times smaller than our resolution. Assuming large optical 
depth, the contribution of this region to the total 7 mm flux density becomes smaller than 
10% of the total flux only if its size is smaller than 2.2 mas (2.6 x 10~^ pc). 

Taking into account these results, we first computed a lower limit to the K band absorp- 
tion from the expected and measured 7 mm flux density, assuming optically thin emission 
for both sources obtaining Ak ~ 1.38 magnitudes, which corresponds to Ay ~ 14 mag. 
An upper limit for the absorption was estimated assuming that the compact source is opti- 
cally thick at 7 mm and does not contribute to the observed flux, resulting in an absorption 
Ay ~ 20 mag. The mean extinction of the optical nebula, derived by Noumaru & Ogura 
(1993) from the Hq;/H/3 ratio is 9.8 mag, implying that there is a strong extinction gradient 
in the north-south direction. 



From the observed size of the extended source, the upper limit for the size of the 
compact component, the upper and lower limits for the absorption, and assuming spherical 



symmetry, we estimated volume densities of (7.4 — 9.8) x 10^ cm~^ and (1.4 — 1.9) x 10^ 
cm~^, representative of compact (CHII) and ultracompact (UCHII) HII regions (Churchwell 
2002), respectively. There is a good agreement between the volume density-size relation (rig 
versus d) for the extended source and what was found for compact and ultra-compact Hll 
regions (Kim & Koo 2001; Martin-Hernandez, van der Hust & Tielens 2003), as can be 
seen in Figure 14. Since we have only an upper limit for the size of the compact source, 
we show in the Figure the track representing the combinations of and d that reproduce 
the observed emission measure. All the points in the track fall inside the dispersion region 
defined in Kim & Koo (2001), which is indicated by the broken hues. 

Finally, we estimated upper and lower limits to the number of Lyman continuum photons 
{N^y) available in both compact and extended HII regions, from the corrected Br7 flux 
densities, using the expression derived by Ho et al. (1990): 



We obtained 2.2 x 10^^ < Niy < 3.9 x 10^^ and 1.6 x 10^^ < < 2.8 x 10^^ photons 
s~^ for the extended and compact regions, respectively. The equivalent single star spectral 
types, obtained from Hanson et al. (1997) are 08.5V-07.5V, and 09V-08.5V, respectively. 
The spectral type of the star associated to the compact region agrees with that calculated in 
section 3.2 from the integrated spectral energy distribution of the IRAS source. However, the 
spectral type of the star expected to ionize the extended region is much earlier than what 
we found for the embedded stars in Figure 3, suggesting that IRSl is the single ionizing 
source of the whole region. Similar results were found by Kim & Koo (1991) for other 
ultracompact cores embedded in more extended HII regions, which they suggested were due 
to champagne flow (Tenorio-Tagle 1979). From the shape of the observed NIR nebula and 
the relative position of the ionizing star (Figure 11), we suggest that this can be also the 



We found that the ionized region classified as the planetary nebula W16-185 is in reality 
powered by a cluster of young and massive stars, still embedded in the parental molecular 
cloud. From the NIR photometry we identified twenty cluster member candidates; the 
majority of them presenting large NIR excess, characteristic of very young objects. 

Using Spitzer infrared images we found that all of selected sources presenting large 





case of W16-185. 



4. 



Conclusions 
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H — K colors have bright Spitzer counterparts. Prom those images we noted the presence of 
at least six highly absorbed mid infrared sources, southeast of the IRSl source. Only one of 
them was detected in the K band of our survey, implying that the extinction must be very 
high in that direction. 

The IRSl source is the most luminous object in the region and also presents large 
infrared excess. This source is the NIR counterpart of the MSX source G326. 7238+00. 6148, 
which is the MIR counterpart of the IRAS 15411-5352 source. From its near to far spectral 
energy distribution and assuming that all the stellar luminosity is used to heat the dust, 
we determined a lower limit for its bolometric luminosity of LboI > 7.7 x 1O*^L0, which 
corresponds to an 08.5V ZAMS star. 

We found very strong Br7 emission at the IRAS source position and more diluted and 
extended emission across the continuum nebula. The first one is very compact and was 
not resolved by our observations, which corresponds to a diameter of less than 0.012 pc 
at the distance of 2.4 kpc, size characteristic of ultracompact HII regions. The size of the 
extended region is 10", corresponding to a linear size of 0.11 pc, typical of compact HII 
regions (Churchwell 2002). 

From the observed Br7 flux we computed the emission measure E and the expected 
free-free emission at radio wavelengths. From the measured 7 mm flux density, we found 
lower and upper limits for the extinction Ay of 14 and 20 mag, respectively. 

From the absorption corrected emission measures, we obtained limits to the volume 
densities of (7.4 — 9.8) x 10^ and (1.4 — 1.9) x 10^ cm~^ for the extended and compact 
regions, respectively. These density values agree with the assumption of the compact and 
ultracompact nature of the sources, as inferred earlier from their sizes. Moreover, the relation 
between volume density and size agrees very well with what is found for other compact regions 
(eg. Kim & Koo, 2001). 

We also obtained upper and lower limits for the number of Lyman continuum photons 
necessary to produce the Br7 emission in both the compact and extended region, obtaining 
2.2 X 10^« < NLy < 3.9 x 10^^ and 1.6 x 10^^ < < 2.8 x 10^^ photons s-^ respectively. 
The spectral type of the single ionizing star in the compact region agrees with that calculated 
from the near to far infrared spectral energy distribution. However, the stars embedded in 
the extended region do not provide enough UV photons to ionize it, suggesting that IRSl 
is probably the main ionizing source of the whole region. In fact, from the shape of the 
observed NIR nebula and from the relative position of IRSl source, we suggest that it could 
be produced by the champagne flow as was found for other ultracompact cores. 
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Fig. 1. — Histograms of J (left), H (center), and K (right) magnitude counts. The complete- 
ness limits values are indicated by the vertical dashed lines. 




Fig. 2. — The K2MASS- Kljsia versus {H — K)l]sia diagram. The dashed hne represents the 
hnear fit to the data, which is shown on the top. We can see that the color term correction 
is very smalL 
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Fig. 3. — The LNA versus 2MASS comparative diagrams. We compared 101 (J), 131 (H) 
and 126 (K) stars in both surveys in an area of about 13 square arcmin around IRAS15411- 
5352. The straight hnes represent the expected relation if the two photometric system were 
equal. 




Fig. 4. — Combined false- image (3.9 x 3.5 square arcmin) made from the J (blue), H (green) 
and nbi^ (red) LNA images. North is to the top east to the left. At the botton-left side 
we present a detailed view (about 1.2 x 1.0 square arcmin) of the near infrared nebula. 
The 5 GHz contour map (Goss & Shaver, 1970) of the region, is shown in the right side 
of the Figure, where the location of IRAS 15408-5356 and IRAS 15411-5352 sources are 
indicated. The IRAS 15411-5352 source is associated with a small deformation in the 5 GHz 
contour map. 
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Fig. 5. — Diagram of spatial distribution of sources detected in the LNA K band image 
(IRAS 1541 1-5352) and for the "control" region taken from the 2MASS survey. The Nebula 
region is indicated by the small box inside the IRAS region. 




Fig. 6. — The (J — H) versus {H — K) comparative diagrams. The locus of the main sequence 
and giant branch taken from Koornneef (1983), are shown by the continuous hues, while the 
two parallel lines (dashed and dotted) follow the reddening vectors given by Fitzpatrick 
(1999). The crosses represent the location of Ay = 0, 10, 20, 30 and 40 magnitudes of visual 
extinction, which are indicated by bold numbers. The dotted line separate the NIR sources 
that present excess emission from that do not. 
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Fig. 7. — The N{H — K) comparative diagram for the control (dotted hne) and Nebula 
(dashed hne) regions. The counts from the Control and Nebula region were normahzed to 
the same area. 




Fig. 8. — The (J — H) versus {H — K) diagram for the nebula region. The locus of the 
main sequence and giant branch taken from Koornneef (1983), are shown by the continuous 
lines, while the two parallel lines (dashed and dotted) follow the reddening vectors given 
by Fitzpatrick (1999). The crosses represent the location of Ay = 0, 10, 20, 30 and 40 
magnitudes of visual extinction, which are indicated by bold numbers. The cluster member 
candidates were selected from the criteria described in the text. The dotted hne separate the 
NIR sources that present excess emission from that do not. Some sources detected only in the 
H and K bands were plotted using the J band completeness magnitude limit (represented 
by arrows). 




Fig. 9. — The 1.4 x 1.0 arcmin images taken from the Spitzer data archive (North is to 
the top, east to the left) of the region around IRAS1541 1-5352. The Infrared Array Camera 
(IRAC) has four infrared bands centered at 3.6, 4.5, 5.8 and 8.0/im with image plate scale 
of 1".2. Some prominent LNA sources are indicated by labels. The small dashed box in 
the southeast of the images delineate a group of highly absorbed mid-infrared sources. The 
small black circle indicate the only one in this group detected in our K band image. 
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Fig. 10. — The Spectral Energy Distribution (between 2.2 and 8.0/im) of the IRS8 (filled 
squares) and IRS 19 sources (filled triangles), together with the SED of a late- type star (open 
circles). The 2.2yum flux densities were taken from the LNA survey; the 3.6-8.0yLtm data were 
taken from IRAC. 
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Fig. 11. — The J versus (J — H) color-magnitude diagram of the soTirccs in Table 1. The 
locus of the main sequence at 2.4 kpc (Giveon et al. 2002) is shown by the continuous line. 
The intrinsic colors were taken from Koornneef (1983) while the absolute J magnitudes were 
calculated from the absolute visual luminosity of ZAMS taken from Hanson et al. (1997). 
The reddening vector for a BOV star (dotted line) was taken from Fitzpatrick (1999). The 
bold numbers indicate the location of Ay — 10, 20 and 30 magnitudes of visual extinction. 
Open up triangles represent the sources that show excess emission in the NIR and filled 
down triangles those that do not. 
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Fig. 12. — K band contour map of the infrared nebula. The contours start at 3.5 x 10~^ 
Jy arcsec"^ and are spaced by 2.1 x 10^^ Jy arcsec"^. The positions (points) and labels 
(numbers) of the selected infrared sources referred in Table 1 are indicated as well as the 
position of the IRAS coordinate error eUipse (dashed line) . 
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Fig. 13. — The spectral energy distribution of the IRSl infrared source. The flux for the three 
near infrared bands (filled circles) were taken from our work. The mid infrared data (filled 
squares) were taken from MSX catalogue (bands A=8.28/im, C=12.13/um, D=14.65/xm and 
E=21.34//m) while the far infrared data (filled down triangles) were taken from the IRAS 
point source catalogue. 




Fig. 14. — The LNA K band image (about 0.45 x 0.5 square arcmin) of the studied region, 
superimposed by the Br7 (white hues) and Ha (black hnes) (given by Noumaru & Ogura, 
1993), contour diagrams. The Br7 contours are flux cahbrated, with the values starting at 
1 X 10~^ Jy arcsec"^ and spaced by 7.5 x 10~^Jy arcsec"^. The two dimensional intensity 
distribution of the Br7 point source, is presented at the right side of the Figure. 
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Fig. 15. — The observed relation between volume electron density Hg and size d for the 
compact (filled rectangle) and ultracompact (continuous line) HII regions. The continuous 
black bold line represents the positions of all points that give the observed emission measure 
for the unresolved Br'j source. The traced line represents the mean relation obtained by 
Kim & Koo (2001) for other compact Hll regions, and the dotted line, the dispersion of the 
observations from which the mean relation was obtained. 
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Table 1. List of the selected near- infrared sources 



IRS 


a(J2000) 


<5(J2000) 


J 


H 


K 





1 


15 


44 


59.46 


-54 


02 


17.4 


11.99 


± 


0.02 


9.77 


± 


0.03 


8.07 


± 


0.03 


2 


15 


45 


00.60 


-54 


02 


05.6 


13.02 


± 


0.03 


11.78 


± 


0.03 


10.71 


± 


0.03 


3 


15 


44 


59.42 


-54 


02 


18.6 








14.11 


± 


0.12 


11.14 


± 


0.10 


4 


15 


44 


57.08 


-54 


02 


30.4 


14.94 


± 


0.02 


13.19 


± 


0.03 


12.15 


± 


0.03 


5 


15 


44 


59.76 


-54 


02 


18.9 


16.76 


± 


0.07 


14.02 


± 


0.08 


12.26 


± 


0.07 


6 


15 


44 


59.90 


-54 


02 


17.2 


16.97 


± 


0.05 


14.49 


± 


0.05 


12.34 


± 


0.09 


7 


15 


44 


58.35 


-55 


02 


33.3 








14.26 


± 


0.05 


12.52 


± 


0.04 


8 


15 


44 


57.14 


-54 


02 


07.5 








17.38 


± 


0.08 


13.00 


± 


0.06 


9 


15 


44 


59.95 


-54 


02 


04.9 


17.03 


± 


0.07 


15.19 


± 


0.04 


13.12 


± 


0.07 


10 


15 


44 


59.69 


-54 


02 


19.9 


16.02 


± 


0.07 


14.46 


± 


0.08 


13.36 


± 


0.06 


11 


15 


44 


59.45 


-54 


02 


09.6 


16.84 


± 


0.05 


14.96 


± 


0.04 


13.53 




0.05 


12 


15 


44 


59.14 


-54 


02 


10.0 








15.59 


± 


0.05 


13.70 


± 


0.06 


13 


15 


45 


00.48 


-54 


02 


02.4 


15.77 


± 


0.05 


14.77 


± 


0.06 


13.80 


± 


0.06 


14 


15 


44 


59.10 


-54 


02 


16.1 








15.87 


± 


0.09 


14.11 


± 


0.06 


15 


15 


44 


58.76 


-54 


02 


10.1 


15.60 


± 


0.06 


14.68 


± 


0.04 


14.12 


± 


0.05 


16 


15 


44 


58.28 


-54 


02 


33.8 


17.57 


± 


0.09 


16.46 




0.13 


14.27 


± 


0.06 


17 


15 


44 


58.47 


-54 


02 


05.3 


16.67 


± 


0.06 


15.79 


± 


0.06 


14.33 


± 


0.06 


18 


15 


45 


01.12 


-54 


02 


02.3 


17.37 


± 


0.09 


16.14 




0.06 


14.44 


± 


0.05 


19 


15 


45 


02.90 


-54 


02 


07.5 








17.79 


± 


0.09 


14.45 


± 


0.08 


20 


15 


44 


59.01 


-54 


02 


12.2 


16.75 


± 


0.09 


15.99 




0.07 


15.33 




0.12 
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Table 2: List of derived parameters 


Parameter 


extended 


compact 


total 


A 1 1 • 1 /51 \ 

Angular diameter ( ) 


10 


< 1.9 




Radius (pc) 


0.06 


< 0.01 




SBr7(obs)" 


12.6 ± 1.6 


9.1 ± 1.2 


21.7±2.0 








5.2 ±0.9 


S43G//2(pred) 


0.85 


0.61 


1.46 


Sj(obs) 






0.17±0.02 


S J (free-free) 






0.10 


Sj^(obs) 






0.29±0.05 


(free- free) 






0.14 


Sx(obs) 






0.38±0.06 


Sa- (free- free) 






0.18 


r(43) 


> 0.003 


> 0.02 






> 5.3 X 10^4 


> 3.8 X 10^6 






7.4 - 9.8 xlO^ 


1.4 - 1.9 xlO^ 




Av 






14-20 




2.2 - 3.9 


1.6-2.8 


3.8-6.7 



"^the Br7 fluxes are in units of 10 ^^erg s ^ 
Hhe remaining fluxes are in Jy 
'^in units of lO^^s"^ 
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